(N 



Probing Planckian physics: resonant production of particles 
during inflation and features in the primordial power spectrum 

Daniel J. H. Chung* 
Randall Physics Laboratory, University of Michigan, Ann Arbor, Michigan 48109-1120 

Edward W. Kolbt 

NASA/Fermilab Astrophysics Center, Fermi National Accelerator Laboratory, Batavia, Illinois 60510-0500, 

and Department of Astronomy and Astrophysics, Enrico Fermi Institute, 

The University of Chicago, Chicago, Illinois 60637-1433 

Antonio Riotto-I- 
Theory Division, CERN, CH-1211 Geneva 23, Switzerland 

'o\' 

G^ ' Igor I. Tkachev^ 

0> ' Theory Division, CERN, CII-1211 Geneva 23, Switzerland, and 

Institute for Nuclear Research, Russian Academy of Sciences, Moscow 117312, Russia 

(September 1999) 

IJ ' The phenomenon of resonant production of particles after inflation has received much attention in 

the past few years. In a new application of resonant production of particles, we consider the effect of 
a resonance during inflation. We show that if the inflaton is coupled to a massive particle, resonant 
production of the particle during inflation modifies the evolution of the inflaton, and may leave an 
imprint in the form of sharp features in the primordial power spectrum. Precision measurements 
►-^ ' of microwave background anisotropies and large-scale structure surveys could be sensitive to the 

r~^ ' features, and probe the spectrum of particles as massive as the Planck scale. 
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0\ ■ I. INTRODUCTION AND MOTIVATION 

CIhI Analysis of the statistics of temperature fluctuations of the cosmic background radiation (CBR) and the distribution 

Jv of galaxies and galaxy clusters in large-scale structure (LSS) surveys yield the presently observed power spectrum of 

fli I scalar density perturbations. The power spectrum is a most useful and fundamental tool for understanding the origin 
r^ t and development of the large-scale structure of the universe. 

It is likely that the observed power spectrum had a primordial origin, resulting from events that occurred in the early 

universe. The primordial power spectrum is not directly observable since it is adulterated by a variety of cosmological 

and astrophysical processes associated with the transition from a radiation-dominated to a matter-dominated universe, 

J-j ' with the decoupling of matter and radiation, with the free streaming of dark matter, and so on. Nevertheless, present 

observations of the power spectrum allow us to infer the form of the primordial power spectrum. 

To a first approximation, present observations of the power spectrum are consistent with it having originated from 
a primordial power spectrum that is featureless, i.e., a primordial spectrum of the form /c" where k is the wavenumber 
and n is the spectral index. However, there are hints of features in the power spectrum from analysis of the Automatic 
Plate Measurement ( APM) survey [|l[ , from apparent regularity in the redshift distribution of galaxies seen in deep 
pencil-beam surveys Q , cluster redshift surveys ||^ , and galaxy surveys Q,^ . 

There are (at least) three possible responses to these "hints." The first response is that the apparent features are 
an artifact of small data sets or data sets with systematic errors, and the effects will eventually disappear with better 
or more complete observations. Another possibility is that the features are real, and they refiect the shape of the 
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primordial power spectrum. A third response is that the primordial power spectrum is smooth, but astrophysical 
processing of the primordial power spectrum produces the features. 

A well known example of the last phenomenon is the "Zel'dovich peaks" in the angular power spectrum of CBR 
temperature fluctuations. The peaks are not believed to reflect features in the primordial power spectrum, but pre- 
sumably result from astrophysical processing, in this case, acoustic oscillations of the baryon-photon fluid. However, 
the features in the matter power spectrum around 100/i^^Mpc are not well fit by acoustic oscillations |g]. 

The leading theory for the origin of the primordial power spectrum is quantum fluctuations during inflation. Simple 
models of inflation produce a featureless, nearly exact power-law spectrum of primordial perturbations. Of course, 
an important qualifier in the previous sentence is "simple." A variety of inflation models that produce features in the 
primordial spectrum have been studied. These include variants of extended inflation [0 that produce rare large voids 
Ipl , models with multiple episodes of inflation 0] , and phase transitions during the inflationary phase 11^] . 

In this paper we propose a mechanism to produce features in the primordial power spectrum by a very simple variant 
of the simplest inflation model. Our basic idea is to couple the inflaton to a massive particle, with resonant production 
of the particle during inflation. Resonant extraction of even a small fraction of the energy of the inflaton field during 
inflation can alter the classical motion of the inflaton and produce features in the primordial power spectrum. 

The primordial power spectrum is most conveniently expressed in terms of the amplitude of the density perturbation 
when it crosses the Hubble radius after inflation, Snik). For reference, the Harrison-Zel'dovich spectrum has Snik) 
independent of /c. A qualitative feel for our mechanism can be obtained by recalling the simple approximation for the 
density perturbations produced during inflation |ll[| , 

6Hik) ^K. (I) 

OTTCp 

Here, H is the expansion rate of the universe and (f> is the velocity of the inflaton field when the comoving wavenumber 
k crossed the Hubble radius during inflation. We find that resonant extraction of inflaton field energy decreases (p, 
leading to an increase in Snik). Because of the resonant nature of the process, the produced feature is rather sharp, 
extending less than a decade in wavenumber. 

There are two reasons why our study may be important. The first reason is that features in the observed power 
spectrum have important astrophysical implications. Fitting the CBR Zel'dovich peaks is an important part of 
cosmological parameter estimation [O]. Spurious peaks in the spectrum arising from resonant particle production 
will complicate the extraction of parameters (although joint parameter estimation from CBR and LSS |1^ may help 
alleviate the problem). 

If baryons are a relatively large fraction of the matter density, the effect of acoustic oscillations on the matter power 
spectrum [|l^,|Ijl may be seen in accurate determinations of the power spectrum in the next generation of large-scale 
structure surveys such as the 2-degree Field project m^ and the Sloan Digital Sky Survey (SDSS) [M, and may be 



part of the parameter estimation program |17|. Again, this program is complicated if the primordial perturbation 
spectrum does not have a smooth power law form p8| . Finally, the program of reconstruction of the inflaton potential 
from CBR and large-scale structure determinations of the power spectrum [|l9|-pi[ may be compromised if there are 
sharp features. 

The second reason we believe our study is relevant is that peaks in the primordial power spectrum may serve as a 
probe of energy scales as large as the Planck scale. We will discuss the particle physics motivation for the coupling 
and masses of the particles in a later section. But to see the possibility of probing Planck-scale physics, simply assume 
a chaotic inflation model with a quadratic potential. The inflaton mass in this model is around 10^^ GeV, which is 
also approximately the expansion rate of the universe during inflation. During inflation the value of the inflaton 
field is a few times the Planck mass, mpi. Since the energy density during the inflationary stage is of the order of 
mim^pi <C TTi^i, in the spirit of effective field theories we may regard the inflaton potential as the low-energy limit 
of some fundamental Planck-scale theory. Furthermore, values of the the inflaton field </> ~ mpi do not necessarily 
lead to a breakdown of the field theory expansion; for example, moduli fields in superstring theories whose vacuum 
expectation values parameterize different string vacua usually have field values of the order of the Planck scale. 

The inflaton may be coupled to other degrees of freedom with mass a few times the Planck mass, but these degrees 
of freedom are usually integrated out and do not appear in the low-energy action. Of relevance for our considerations 
is the possibility that the inflaton field couples to a particle of mass a few times the Planck mass. Normally, the 
massive particle would be integrated out and not affect the low-energy theory. In the inflaton models we study, 
observable length scales exit the Hubble radius during inflation when the value of the inflaton field is a few times the 
Planck mass. If the inflaton couples to a fermion of mass m (assumed to be of order of mpi) with Yukawa coupling 



A, then the ^-dependent mass of the fermion is M{(j)) = m — A(/>.[| The cj) ground state is (/> = 0, and in the ground 
state the fermion is supermassive and not part of the effective low-energy theory. But when cj) ~ m/A, the effective 
mass of the fermion vanishes, and it cannot be integrated out. When ~ (/>» = rn/X and the fermions are hght, and 
they can be efficiently produced |22f| . Since ip > mpi during inflation, it can act as a probe of Planckian physics. 

Ours is not the first study of parametric creation of fermions. Pure gravitational creation was considered in Refs. 
[^,|^. Creation by an oscillating background field in Minkowski space was considered in Refs. |^,Q. Resonant 
production of massless fermions after inflation (i.e., during preheating) in a X(fA inflation model was studied by 
Greene and Kofman |g^ . Non-thermal productions of gravitinos have been considered in Ref . |2j . Finally, resonant 
production of massive fermions during preheating a/ier inflation in an m'^(f>'^ inflation model was considered by Giudice, 
Peloso, Riotto, and Tkachev |g^. We draw heavily from the results in these important papers. The new twist we 
present is a change in the venue for particle production: we assume the resonance occurs during inflation. We believe 
this may be the most easily tested possibility because it leads to directly observable effects in the power spectrum. 

In the next section we will provide more details about the couplings in the model we study and discuss the basics of 
resonant particle creation during inflation. In Sec. Ill we consider resonant production where back reactions are not 
important. We make some analytic approximations leading to an estimate of the effect on the density perturbation 
spectrum. In Sec. IV we present the results of a numerical calculation of the perturbation spectrum including the 
effect of resonant particle creation in the regime where back reactions are important. In Sec. V we discuss the prospect 
of detection of the type of features we produce. Finally we present a brief summary in the concluding section. 

II. HEAVY PARTICLE PRODUCTION DURING INFLATION 

Resonant particle production can occur for a range of choices of inflation models as well as different heavy particle 
masses and couplings of the massive state to the inflaton. We will develop analytic techniques and present numerical 
results for one model, and mention the generalization of our results where appropriate. 

We will assume that the inflaton potential during inflation is given by a simple quadratic potential of the form 

F(0) = ^ml4>^ . (2) 

The value of the inflaton mass can be fixed by normalization of the perturbation spectrum, with result to^ = 10^^ GeV. 
Inflation occurs for Planckian values of the inflaton field, cj) Z mpi, where mpi is the Planck mass. Inflation ends 
when 4) ~ 0.2topl, and scales relevant for CBR fluctuations and large-scale structure studies cross the Hubble radius 
during inflation when is in the range of approximately 2mpi to Smp; . 

Now we will couple a fermion ij) of mass M to the inflaton by a Yukawa term of the formQ 

Cy = -XHil^ . (3) 

If (j) is nonzero, then the effective mass of the fermion i^ 

M{4>) ^m~ \(t) . (4) 

These choices lead to a critical value of the inflaton field, which we will denote as 0,, where the effective mass of the 
fermion vanishes: 

(/)» — to/A . (5) 

Resonant fermion production will occur in a narrow interval around (f) — (f>^, 123] . 



'^The important feature here is that the 0-dependent mass vanish during inflation. This can be accomplished by choosing a 
negative sign for the Yukawa term, or by taking to be negative. In this paper we will choose positive value of (f> and a negative 
Yukawa term. An important notational convention is the distinction between m and M[(l)) — m ~ \(j). 

^Here we treat a Dirac fermion, but there is no qualitative difi'erence if one chooses Majorana fermions. 

^Notice that radiative corrections from the heavy fermion may spoil the flatness of the inflaton potential. As usual, one can 
invoke supersymmetry to preserve the flatness. If this is the case and the starting superpotential is 214^ = m^(lP' + {m — \(f>) ip^ , 
it is easy to show that radiative corrections are negligible, with the one-loop effective potential S2tt^Vi ~ X^m^tj)^ log(A^cjf>^/A^), 
where A is the renormalization scale. 



Although we chose a fermion for the massive state, one might imagine that the massive state is a boson. Recall 
that the requirement for particle production is that the effective mass term goes through zero at 0* . If the effective 
mass term goes through zero, then one might expect that it is negative either above or below 0*. Since the sign of 
a fermion mass term is irrelevant, a negative mass term is not a problem. If the massive state is a boson, then the 
quadratic nature of the mass term implies one has to worry about tachyonic modes. So if the massive state is a boson, 
then the effective mass term should remain positive both above and below (jj^. This can be done, for instance, by 
taking a boson a with potential V{a) — a^ {rria — gcj)) + Ao-ct'*, where Too- > 0. The vacuum expectation value of the 
field a vanishes. The effective mass term of the boson, Mj = 2 [rria — g(f>) , remains positive definite and vanishes at 
(f>* = rricr/g. 

Now consider fermion production. In an expanding universe, the fermion field ^ satisfies the Dirac equation in 
conformal time rj (adrj = dt) , 

-^ + 4^7°-M| ^^ = 0, (6) 

a 2 J 

where a is the scale factor, 7'^ are the flat-space gamma matrices, the spatial derivatives are with regard to comoving 
coordinates, and the Hubble expansion rate is H = a/ a = a' /a? (here and throughout the paper a prime superscript 
implies d/drj) . This equation can be reduced to a more familiar form of the Dirac equation by defining a new field 
X = a~^'^4'- In terms of the new field Xj the Dirac equation in an expanding universe becomes 

[*^-a(77)A/(77)]x-0. (7) 

Of course M depends on 77 through the dependence of (j) on rj. 

The field x can be expanded in terms of Fourier modes of the form 

^(^) = / TTl^ ^'""^ E [Mk,v)ar{k)+v,.{k,r^)bl{-k)] , (8) 
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where the summation is over spin and Vr{k) = Cu^{—k). The canonical anticommutation relations imposed upon the 
creation and annihilation operators may be used to normalize the spinors u and v. 

Defining Ur = [u+{ii)ipr{k),ru-{ri)4'r{k)] and Vr = [rv+{ri)-4!r{k),V-{rii)'ipr{k)] , where iprik) are the two- 
component eigenvectors of the helicity operators, and using a representation where 7° = diag(]l, — 1), Eq. (^) can 
be written as two uncoupled second-order differential equations for u+ and u-: 

< + [luI ± i{aMy] u± = , (9) 

where, ^^(77) — k'^ + AP {r])a'^ {rj) and M{r]) — m — X(j){-q). Notice that the normahzation conditions ulug — v^Vg — 2Srs 
and ulvg = are preserved during the evolution. Furthermore, sometimes it is useful to remember that w+ has the 
same evolution as — ul and w_ has the same evolution as u'^. 

In order to calculate the number density, we must first diagonalize the Hamiltonian. In the basis of Eq. dq) the 
Hamiltonian is 



H{rj) = / d'kJ2{Ek{v) K(fc)«r(fc) - br{k)bl{k)] + Fk{v)br{-k)ar{k) + F*{7j)al{k)bl{~k)} , (10) 

•' r 

where the equations of motion can be used to express Ek and Fk in terms of u^ and W-:Q 



Ek = kRe{u*^u^) + aM ( 1 - |u+ 

k 
Fk ^ -{u\- ui) + aMu+u^ . (11) 

In order to calculate particle production one wants to write the Hamiltonian in terms of creation and annihilation 
operators that are diagonal. To do this one defines a new set of creation and annihilation operators, d and b\ related 



*Here we choose the momentum k along the third axis and use the representation in which 7'^ = ( 



to the original creation and annihilation operators a and b^ through the (time-dependent) Bogolyubov coefficients ak 
and pk, 

a{k) ^ ak{v)a{k) + PkivPi-k) , 
b^{k)^-f3Uv)aik)+aU7j)b\-k). (12) 

The Bogolyubov coefficients will be chosen to diagonalize the Hamiltonian. Using the fact that the canonical com- 
mutation relations imply jafcp -I- |/3fcp = 1, the choice 

ak _ Ek+ ujk |a |2 _ |JfcP ..„^ 

Pk^ F^ ' 1^'=' -2u;kiu;k+Ek) ' ^'^^ 

results in a diagonal Hamiltonian, 

H{t^)= f d\J2LUkiv)[alik)arik) + bl{k)brik)] . (14) 

We define the initial vacuum state |0) such that a|0) = 6|0) = 0. The initial conditions corresponding to the 
no-particle state are 



, , ojk T Ma , , s , ^ , s , s 

u±{0)^\-^— ; w±(0) = ifcu^(0)T«aMu±(0) . (15) 

V ^k 

The (quasi) particle number operator can be defined as JV = al{k)dr{k), and the particle number density n is 
(including the two degrees of freedom from the spin) 

1 f^ 

n{rj)^{0\M/V\0) = ^^-- dkk'\(3k\^. (16) 

TT-^a-^W Jo 

An equal amount of antiparticles is produced. 

Including the coupling of the fermion to the inflaton, the inflaton equation of motion in the Hartree approximation 

is 

4> + 3H^+^-NX{^i;) =0 , (17) 

dip 

where we have generalized to the case of N fermions of bare mass m coupled to the inflaton. 

The product (ipi^) can be expressed as momentum integration of the mode functions using the field decomposition, 
Eq. (H). A straightforward averaging leads to ultraviolet divergences, therefore, the quantity (ipip) must be regular- 
ized. As in the case of Minkowski spacetime, the regularization amounts to normal ordering or, equivalcntly, to the 
subtraction of zero-point vacuum fluctuations. To obtain a finite result, it is necessary to express the operator ^p'lp 
in normal-ordered form and subtract the part due to vacuum fluctuations. The normal-ordered V'V' operator has the 
form 

N^ {^i:) = ^^ - (O^li^^PlO^) , (18) 

where the vacuum |0^) is defined by 

a|0„)=6|0„)=0. (19) 

The vacuum averaging in Eq. (O) is defined as averaging with respect to the original vacuum state (we remind the 
reader that we are working in the Heisenberg representation) 

(V^V) = (0|iV, iiJ^) |0) . (20) 



Using the operator redefinitions given in Eq. (12), it is straightforward to show that 



W) ^J^Jd'k [(|«+P- 1) \l3kf+Re{ak(3luXu*_)] . (21) 

Making use of the formulae in Eq. dlq), we arrive at the following expression for the regularized average of -tpip: 



<^"^^ = 7^ fd'k(\u^\' + ^~l]. (22) 

[2'Kay J \ ujk ' 



In the case of Majorana spinors the numerical factor is half that of Eq. (ES 

We will use Eq. (|2^) when integrating the equation of motion [given in Eq. (|l^)] for the inflaton field. We have 
also consistently included the contribution from (V'V') in the equation of state when integrating Einstein's equations 
for the scale factor. We integrated the system of all equations in the conformal frame with Eq. (||) the equation of 
motion for the mode functions. For the case of an inflaton potential given by Eq. (ph, Eq. (p7|) takes the form 



-,2 2 ^ 



if" + I m^a^ 1 (^ - a-^NX{iji}) = , (23) 

where f = acj). 

At this point it is useful to discuss the unknown parameters of the calculation. First, we do not know the inflaton 
potential {V). We simply assume it is a quadratic large-field model for the purposes of detailed calculations and 
return to this issue in a later section. We do not know the mass of the fermion {m) and the Yukawa coupling of the 
fermion to the inflaton (A). We will assume that m and A are such that i^*, defined in Eq. (g), falls in the range of 
(j) that corresponds to the epoch of inflation when astrophysically accessible scales are crossing the Hubble radius. In 
the large-field inflation model we study, this corresponds to m/A a few times the Planck mass. Particle production is 
most efficient when A is large, so m also must be large to have (/)* in the relevant region. Finally, we do not know the 
number of fermion fields that couple to the inflaton {N). This number is potentially very large because there are often 
fermion representations of large dimension in superunified theories. So with the assumption of the large-field quadratic 
inflation model, the parameter space is spanned by {m, A, A^}, with the ratio w/A reasonably well constrained. 

III. ANALYTIC APPROXIMATIONS AND NUMERICAL RESULTS WHEN BACK REACTIONS ARE 

UNIMPORTANT 

If only a small fraction of inflaton energy is extracted in resonant production of particles, then it is possible to find 
reasonably accurate analytic approximations for the features produced in the primordial perturbation spectrum. 

There are three steps in obtaining analytic approximations for the modification to the perturbation spectrum due 
to resonant particle production during inflation. The first step is to determine the efficiency of resonant fermion 
production. The next step is to estimate the effect of fermion production on cj). Finally, it is necessary to convert a 
change in (j) to the modification of the perturbation spectrum. 

A. Fermion Production 

First consider the question of particle production. Start with the wave equation, Eq. (0): 

u± + {k^ + M'^a^)u± ± i{aM)'u± = 0. (24) 

Now define the (conformal) time 77* to be the time when M{ri^,) — 0. Recalling that M = m — X(j), in the vicinity of 
77 = 77* we may expand Mijf) as 

Af(77) = A0:(77-7y,) + .-- , (25) 

where 4>'^ = 4>' {"q^,). Now let's assume that resonant production happens in a very narrow region of conformal time 
and the scale factor does not change much during particle production (i.e., a' — 0). Choosing a* = 0(77*) — 1, then 
around rj^ the wave equation is 

ul + [fc2 + A20t'(7; - 77,)2] u± T i\4>'.u± = . (26) 

Now defining new variables p = k/ ■\/A(/>'^ , r = -^/At/i'^ (77 — 77*), and denoting d/dT by a dot, the wave equation becomes 

u± + {p^Ti + t"^) u± = . (27) 

The solutions to this equation are the parabolic cylinder functions. Choosing the boundary conditions associated with 
the no-particle initial state, one can estimate the phase-space density of the created particles to be (cf. the Bosonic 
case discussed in Refs. p9|-pl|) 



\h 



exp - 






(28) 



The argument of the exponential appearing in Eq. ( p8| ) can be expressed in terms of a physical momentum, fcp, the 
time derivative of the inflaton field, and the expansion rate at resonance, as 






A \6J 






(29) 



Here H^ is the Hubble expansion rate and (/>* is the unperturbed velocity of the inflaton field at 77 = 77, 
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FIG. 1. The spectrum of produced fermions around rj 
analytic approximation. 



rj, as a function of physical momentum kp. Also shown is the 



If the resonance occurs at 



2m PI, then H^ ~ 4toj, and 



Q.lQmpimM. The crucial fact is that in this model 



the ratio Hl/\(f)^\ is very small, about 10^**. This allows the phase-space density to be large in regions where the 
physical momentum is larger than the inverse of the Hubble radius. The smallness of this ratio should obtain in all 
slow- roll inflation models. 

Fig. n^ shows an example of the comparison between the analytic formula of Ec^. (Eq) and a numerical solution of 
the equations. The momentum is in units of kp/H^, where k„is the physical momentum at the time 77,. 

There are two important results one may learn from Fig. nl The first is that the analytic approximation is quite 
accurate. The second important fact is that, as expected from the discussion above, the phase-space density is large 
for k greater than i?*.Q| This means that the effect is dominated by sub-Hubble-radius wavelength modes. 

The physical number density of each Dirac fermion species (particles plus antiparticles for a total of four degrees 
of freedom) produced at 77* is given by 



n{r]^ 



dkpkl\(3k\ 




(30) 



For rj > 1]^,, the number density of fermions simply decreases as a '^. For r] > rj^, the fermion mass is approximately tti, 
so knowing the number density at 77 = 77* allows us to estimate the contribution of the produced fermions for 7; > 77* . 



The contribution per logarithmic interval in k is proportional to k \f3k\ , so the peak contribution to the number density of 
created particles occurs for kp ^ H,. 



B. Inflaton Velocity 

Now we turn to the issue of how resonant production of the fcrniions modifies <j). Let's start with the equation of 
motion for the inflaton field, Eq. (17). Here we are interested in change in the velocity of the inflaton field caused by 
particle creation. Consider a small region around the time of particle creation, t = t^. In this small region H = H^, 
and dV/d(j) = {dV/d(j))^ are approximately constant. 

We may also approximate (V'V') by assuming the number density of fermions is zero before i*, equal to n, = n(?7») 
[given by Eq. (|0)] at t = i*, and decreases as a~^ thereafter: 



{^^P) = n, d{t ~ U) (— ) ~ n, 0{t - U) exp[-3ff*(t - U)] 



(31) 



Here for simplicity we assume an exponential growth of the scale factor during inflation. 

The fact that {"ipip) can be approximated by n* may be understood as follows. The mode functions have an adiabatic 
evolution before and after the time 77* when M{ri^,) = 0. The nonadiabatic changes of the mode functions occur only 
in the vicinity of 77* . Therefore we expect the mode function 



u+(77 < r]^,,k) 



Ma 



exp 



Uk dri 



(32) 



to be the adiabatic solution of Eq. (H) before the particle production at time 77*. After the production time, the 
mode function u+ is again an adiabatic solution, with the difference that it picks up a negative frequency. It may be 
expressed in the form 



u+{ri > r]^,,k) = afe\/l exp 

^k 



-I I LUk df] \ - /5feW 1 H exp 



LUk drj 



(33) 



Inserting this expression into Eq. (22) and making use of the property \ak\ + |/3fe| = 1 , we obtain 



i^iP) 



(27ra)3 



,0 , I Ma , „ ,2 fc „ 



akPl exp -2i I ujk drj 



(34) 



The Bogolyubov coefficient |/3fc| may be well approximated by unity for k ^ k^: — ^J\ \ (j)'^ \ and zero otherwise. 
Furthermore, at 77 > 77* it is safe to neglect the inflaton field dependent part in the mass Af of the fermion 7/; and we 
have M ~ 777 ^ fc*. Under these circumstances, Eq. (34) reduces to 



W) = 



(27ra)3 



d'k\—\l3k\ 

OJk 



(27ra)3 



d^k\(3kr=n,. 



(35) 



This proves that (ip'ip) can be approximated by 77*. 

With this approximation for {-ipip), the first-order differential equation for [ Eq. (|l7|)] is 



dt 



+ 3ff*0 + {dV/d(j))^ - NXru6{t - U) exp[-3iJ*(i - U)] = 



(36) 



The solution of this equation for i > t* is 

^(t > U) = 4>^ exp[-3H4i - t,)] - '■ /„ ''* {1 - exp[-3i?*(i - t,)]} + NXn^{t - U) exp[-3iJ*(< - U)] . (37) 

3/1* 

We may easily calculate the change in (j) due to particle creation: 



A(t>(t > Q ^ (j)(t > U) - (j)(t>U) = NXnJt~U) exp\-3HM-U)] 

I J A=0 



(38) 
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FIG. 2. The numerical result for the evolution of tp, including the effect of resonant particle production, is shown by the 
solid curve. Also shown by the dashed curve (nearly indistinguishable from the numerical result) is the analytic expression of 
Eq. ph. 



(39) 



The maximum size of the feature can also be estimated. Equation ( pSD has a maximum at i = (3i?* ) of 

max[A0] _ A0(t = l/3iJ*) _ NX^/^ |0*|3/2 



mpim^ 



mpim^ 



BcTT^ H^mpim^ 



where we have used Eq. ( pOJ ) for n*. 

In Fig. we present the result of a numerical calculation. Note that m.ax[A(/)] is about O.OSmpiTTT.^. For the model of 
Fig. 0, H„ ~ 4m0 and </)* ~ O.lGmpjm^, and from Eq. (39) we predict max[A0] ~ O.OSmpim^, in excellent agreement 
with the numerical results. Not just the depth of the feature, but the entire shape of the feature is well approximated 
by Eq. (H). 

We have assumed that the fermion is stable, or at least has a lifetime long compared to the time it takes to form 
the feature, t — t* ~ 0.2m7^. This is not unreasonable at all. Indeed, in spite of the fact that the bare mass of 
the fermion is much larger than the inflaton mass, the fermions are produced when they are massless. Suppose we 
parameterize the decay rate of the fermion field by F = aM , where a is a perturbative constant smaller then unity. 
Using the expansion of Eq. (K^) and imposing the condition F > iJ, we estimate that the fermions decay at time to 
such that to — t^ c^ (Xacf)^) . The fcrmionic decays take place after the formation of the features in the spectrum if 
to — t*^ O-'^iTT'd) ■ This amounts to requiring a ^ (5/A)(m0/(/)*), which is not very restrictive. Another option is that 
the perturbative decay rate F identically vanishes for kinematical reasons at the inflationary stage. This will happen 
if the particles coupled to the fermion ip are heavier than the mass of the ^ during inflation. 

In the opposite extreme, F <C -ff*, we may assume that the fermions decay very soon after they are produced. 
In this case, the factor of exp[— 3iJ*(i — t*)] in Eq. ( pi| ) would be modified. Again using the fact that around the 
resonance M ~ \(f>^,{t — t*), the factor becomes exp[— 3iJ,(i — i*)] exp[— aA(/),(t — i*)^]- One can follow through and 
derive a similar analytic approximation to A0 [cf. Eq. (p8|)], 



A(p = NXn^, exp[— 3i?,(t — t,) 



1 



■ erf 



aA0, 



aX(f>t.{t — i,) 



(40) 



Of course in the limit that the particle is stable, a ^ 0, the error function may be expanded to recover Eq. (pS 
The decay of the massive state results in a feature that is narrower and less pronounced than if the massive particle 
energy density decreases simply as matter (see Fig. m. The feature in (p will map onto the feature in Sh as discussed 
in the next section. 

We do not know the relevant value of a, but if we choose a random value, say a — (137)^^, then y aA0, = 34 for 
A = 1. This would make a few percent change in the amplitude of the feature and make it a bit narrower (see Fig. 0). 
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FIG. 3. Analytic solutions for </> assuming no decay of the fermion (the solid curve) and with decay parameters 

\/a\(fit,/m^ = 10^, 10^, and 10^. The feature disappears as yaXcpt/m^ -^ oo, and as the quantity approaches zero, we 
recover the result of no decay. 



C. Effect on the Perturbation Spectrum 



The simple expression for the perturbation amplitude in terms of H and given in the introduction, Eq. p), 
must be modified for several reasons. First of all, a correct treatment of the perturbations produced during inflation 
would involve a variational calculation |32f| . The results of this formalism will be reported in another paper. The 
other reasons the simple expression must be modified can be seen by recalling the derivation of Eq. (11]). Ignoring 
normalization for a moment, the correct starting point is |33[| 



Sh oc 



5p 

P + P 



(41) 



= (tP , but around the 



First consider the denominator of Eq. (fil\). If only the inflaton contributes to p and p, then p+p = 
resonance the fermions contribute to p and p. An expression for p + p that includes this effect is p + p 

In the absence of fermion production, the numerator of Eq. (ETl) is just 6p — 6(p{dV/d(p), with S(p given by quantum 
fluctuations of the nonzero mode of the inflaton field, S(p = HJ2tt. However, fermions produced by the resonance 
contribute to the fluctuations for two reasons: fluctuations in the fermion field lead to fluctuations in the energy 
density, and coupling of the produced fermions to the nonzero modes of the inflaton field affect the perturbations. 
Since the zero mode of the inflaton field still dominates the energy density, these effects are not expected to be large. 
If they are included, they will increase dp, hence increase the size of the feature in the primordial power spectrum. We 
will ignore the extra contributions to 6p here, and simply assume that 5p = 6(j){dV/d(j)) = {H/2TT){dV/d(j)). Finally, 
one should not use the slow-roll equation (dV/dcj) = —3H(p), because may be large when fermions are produced. 

In spite of the above considerations, we find that use of the simple expression of Eq. (|l|) allows reasonably accurate 
description of the feature in the primordial power spectrum. In the discussion below we describe both analytic and 
numerical calculations using the following expressions for Sh'- 



Sh = 



Hi 

57r0 

4:H{dV/d(j)) 

15Hm% 



analytic 



numerical 



(42) 



For the analytic approximations we use H^ = (87r/3)F/mp; and 



determined by Eq. (^ 
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With the approximation that the change in cj) directly translates into a change in the perturbation spectrum, 
the character of the features produced in the spectrum can be understood from Eq. (|38|). Since A(j) decays with a 
characteristic width of 3iJ*, the width of the feature will correspond to less than one e-fold of expansion. For a given 
change in (j) we can estimate the effect on the density perturbation. For the moment, assume the simple result for 
the density perturbations of Eq. p). For reasons discussed above, we expect the width of the feature to be less than 
one e-fold of expansion. Defining a* to be the value of the scale factor at resonance, a numerical calculation of the 
spectrum of density perturbations is shown in Fig. Has a function of a/ a*. As expected, the width of the feature is 
less than a single e-fold, and the peak amplitude of the feature corresponds to a peak change in (j) estimated from Eq. 




FIG. 4. Resonant particle production produces a peak in the perturbation spectrum as shown in the figure. The solid 
curve is the numerical result, the dashed curve is the analytic form of Eq. (ksl). Finally, the dotted curve indicates the power 
spectrum in the absence of resonant particle creation. 



An analytic expression for Snik) can be obtained using the analytic expression for 
mated by a function of the form 



Snik) 



[SHik)] 



x=o 



1 — 6{a — a,)|(/),| ^N\n^,H^ "^(a,/a)3 ln(a/a,) 



The feature can be approxi- 



(43) 



This expression is also shown in Fig. H. Clearly it is a very reasonable approximation. 

In the comparison between analytic and numerical results presented in this section we have assumed a single fermion 
species coupled to the inflaton. One might imaging a large number of fermion species coupled to the inflaton. In grand 
unified theories (GUTs) based upon gauge groups such as S'O(IO), -Eg, and so on, representations of large dimension 
are easily obtained. 

If back reactions are ignored, then the effect of fermion production on the perturbation spectrum is linear in N. 
However, for large values of iVA^/^ the back reactions are important. We turn to such a case in the next section. 

Before concluding this section we restate the important results: The feature in the spectrum should be about one 
e-fold in width, and the amplitude of the feature should scale as NX^/'^ until back reactions become important. 

Finally, we note that we can estimate how the change in the perturbation spectrum scales with the inflaton potential. 
Recalling that 5h = H'^/5'^4>*, then 



SH{k) 



A(j)^ 



NX^/^ 31/4 
6e7r3 83/4 ^, 



1 



3/4 



mli{dV/d(l)), 



V, 



3/2 



1/2 



(44) 



where here we have used the slow-roll equation of motion, 3iJ*^* — 
recall the approximation 



'{dV/d(j))* and used H"^ — 87rK./3mp;. If we 
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OH 



1/2 



y3/2 



'p; 



dV/d<j) 



(45) 



we see that the term in parenthesis in Eq. ( [44| ) is simply proportional to SJj . Therefore, if dV/d(j) and V are adjusted 
for a given 6h, then the relative amplitude of the feature will be independent of V and dV/d4>. This implies that the 
results of the paper can be applied to all single-field inflation models, not just the large-field chaotic inflation model 
we used. 

IV. NUMERICAL RESULTS WHEN BACK REACTIONS ARE IMPORTANT 

We do not expect a qualitative difference in choosing 0* ~ 3mpi , which corresponds to scales relevant for the CBR, 
and a lower value of (/)* . Numerically it is convenient to work with smaller 0* ; in our runs we choose ^* ~ 2mp; . 

Since we already presented results for the case A = 1, m = 2 x lO^^GeV, and A^ = 1, for completeness we also 
present in Fig. ^ the results for (^(i) and H{t) in that model. 

Note that resonant fermion production has only a marginal effect on the evolution of cj) and H, but of course it has 
a large effect on 0. Near resonance, (/)* ~ 2mpi and iJ* ~ ^m^. The resulting feature in the perturbation spectrum 
was shown in the previous section. 




-0.1 0.1 0.2 ^0.3, 0.4 0.5 0.6 0.7 0.8 

(t-tj m^ 



FIG. 5. The evolution of H and (f> in this model is not noticeably altered by resonant fermion production. Only small 
inflections in (f){t) and H{t) are noticeable around the time of resonance, t ~tt where t, is the time when (f) — (f),. 

Now let's study a model with an even larger amount of fermion creation, where back reactions are quite important. 
Let's consider the case of A = 0.316, N = 100, and m — 6.3 x 10^* GeV. As in the previous example, the resonance 
occurs for 0* ~ 2mpi. The evolution of (j> and H in this model is shown in Fig. 0. Clearly the evolution of both 
and H are modified by resonant particle production. The change in (f) is even more dramatic than in the previous 
example. 
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FIG. 6. The evolution of H and in this model is noticeably altered by resonant fermion production. The flattening in the 
evolution of (j) and H begins at f = f, and continues for a time of about 0.2m~^ , after which (f> and H resume decreasing at the 
rate they would if resonant particle creation had not occurred. 



The estimate of the change in in the previous section overestimates the true change. The scaling of Eq. (pE 
predicts max[A0]/mp;r7i0 of 100(0.316)'^/^ times the previous estimate of 0.03. This would predict a change in 
of 0.17, which in fact is larger than (j>^,. The back reaction of the created fermions limits the value of max[A0]. A 
numerical calculation of 4> in this model is presented in Fig. R. While the change in (/) is limited by the back reaction, 
nevertheless it is quite large. 




-0.1 0.1 0.2 ^0.3, 0.4 0.5 0.6 0.7 0.8 

(t-tj m^ 



FIG. 7. The change in (p in this model is modifled by the back reaction limiting the energy extracted from the inflaton fleld. 
Also shown by the dashed line is the analytic expression of Eq. (M), with A fit to agree with max[A(;/)]. 



Even though the value of max [ At/)] is modified by the back reactions, the feature still has the form 

A0(i >t^)=B{t-U) exp[-3iJ,(t - t,)] , 
where _B is a constant, B < NXn^,. The inequality is saturated in the absence of back reactions. 



(46) 



13 



A numerical calculation of the primordial perturbation spectrum is shown in Fig. ra as a function of k/k^, 
wavenumber k is related to a by 
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Oo-ffo 
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end 



v: 



1/4 



end 






A physical 



(47) 



where the subscript '0' indicates present values; the subscript '/c' specifies the value when the wave number k crosses 
the Hubble radius during inflation (fc — aH); the subscript 'end' specifies the value at the end of inflation; and p^ch 
is the energy density of the universe after reheating to the standard hot big bang evolution. The various values of V 
denote the potential energy density at the indicated epoch. This calculation assumes that instantaneous transitions 
occur between regimes, and that during reheating the universe behaves as if matter-dominated. 

Clearly the exact relationship between a/a^, and a presently observed value of wavenumber is model-dependent. 
Although the location of the peak in wavenumber (denoted by /c* ) , cannot be unambiguously stated, the width of the 
feature in wavenumber is proportional to the width of the feature in the ratio a/a^, i.e., fc/fc, oc a/a^. 

Also shown in Fig. g for comparison is a Harrison-Zel'dovich (H-Z) spectrum (spectral index for the primordial 
spectrum of n = 0), and a numerical calculation of the spectrum [assuming the spectrum is given by Eq. (m)] in 
the model without resonant fermion production. The effective spectral index in the nonresonant inflation model is 
ns = l + d\nSjf/d\nk:^0.96. 

This model produces a huge spike in the primordial power spectrum, more than a factor of three times the underlying 
spectrum. The feature is fit by a function of the form of Eq. ( [43[) : 



fe(fc) 



[SH{k)],^, 



A=0 



l-0{k- K)^.lhA{K/kY \n{k/k^ 



(48) 



If back reactions are ignored, then A = |0,|^^A^An*iJ^^/8.15 ~ 0.23A^A^/^. The coefficient A is limited by back 
reactions to be such that the denominator is positive for k > k^,; this implies A < 1. The first model {N = 1, A = 1) 
is unaffected by back reactions, and A = 0.23. Since back reactions are important in the second model {N = 100, 
A = 0.316), the numerical fit (A = 0.7) is smaller than expected from the analytic expression A = 0.23NX^/^ = 1.3. 
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FIG. 8. Resonant particle production produces a peak in the perturbation spectrum as shown by the solid curve in the 
figure. Shown by a dashed curve is the analytic form of Snik) in Eq. (WSl) with the coefficient A chosen to give the best 
fit. Finally, also shown is a Harrison-Zel'dovich spectrum and the spectrum in the inflation model without resonant particle 
creation. 

V. DISCUSSION AND PROSPECTS FOR DETECTION OF FEATURES 

In this paper we have described a reasonably straightforward and simple mechanism to produce features in the 
spectrum. The features can be characterized by a central location in wavenumber, a width, and an amplitude. 
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The location of the peak in wavenumber depends on the value of (/), — m/\. In the large-field inflation model 
we study, (/>* has to be of order two to four times the Planck mass to have any hope of being in the observationally 
accessible range. Since to produce a noticeable peak A cannot be too small, m must be of order mpi. 

The narrow width of the feature, less than about a decade in wavenumber (corresponding to about a single e-fold 
of inflation), is reasonably independent of the model parameters. 

The amplitude of the feature does depend on the model parameters, in particular on the combination NX^'"^ . In 
the first model we presented (A = 1, iV = 1, tti = 2 x 10^^ GeV), the amplitude of the feature was a factor of 1.3 larger 
than the underlying spectrum. The second model (A = 0.316, N = 100, m = 6.3 x 10^^ GeV) had a peak increase 
in the power spectrum of about 3.4 times the non-resonant result, or about 2.6 times the first model. Scaling with 
iVA^/^ would predict a peak increase of the second model of 5.6 (rather than 2.6) times that of the first model. The 
failure of the peak amplitude to scale as predicted is because of the back reaction limiting fermion creation. 

Our conclusion is that the characteristic signature of our mechanism is a narrow spike in the power spectrum of 
width less than a decade in wavenumber and amplitude anywhere from near zero to factors larger than three. An 
exceedingly good analytic form for the shape of the feature is given in Eq. (^). The peak amplification is [1 — A)^^ 
(recall that A <1). The first model had A — 0.23, while the second model had A = 0.7. 

We now turn to a discussion of how such a feature might be detected in near-future experiments that will measure 
the present power spectrum. Since the feature produced by resonant particle production is approximately described by 
a single parameter (0 < A < 1), we can estimate the sensitivity of large-scale structure surveys and CBR experiments 
to a resonant particle production feature. The limit on A of course will depend on the value of fc». 

It is a relatively straightforward exercise to see how a feature in the spectrum would appear in determination of the 
power spectrum from large-scale structure surveys. The first step is to convert the primordial spectrum, (5//(fc), which 
is the amplitude of the perturbation as wavenumber k enters the horizon, to the present-day power spectrum, P{k), 
which describes the amplitude of the perturbation at a fixed time. This is done by means of the transfer function, 

r(fc) @: 

^\p{k)^(^\THk)5l{k). (49) 



27r2 ^ ' \aHa 

The transfer function depends upon a whole range of cosmological parameters (51, Hq, A, Qb, and so on). This 
expression should be valid in the linear regime of the perturbations, which in comoving wavenumber is approxi- 
mately k < 0.2/iMpc~^. Features in Snik) are directly passed through to P(fc), and in the linear regime should be 
straightforward to identify. 

For a volume-limited survey, the uncertainty due to cosmic variance and shot noise in the estimated power per 
mode is f^5|J36| 



AP(fc) ^ / (2^)3 
P{k) V ^sVk 



nP{k) 



(50) 



where n is the mean density of the sample of survey volume Vs and the estimates average over a shell in Fourier space 
of volume Vk = Airk'^Ak. The width of the bins is Afc = 2TrR~^, where it! is the depth of the survey. 

Consider the sensitivity of determinations of the power spectrum expected from the bright red galaxy (BRG) 
subsample of the SDSS. The volume-limited BRG subsample will include approximately 10^ galaxies in volume of 
effective depth Ih^^ Gpc covering tt steradians. It is expected that BRGs are biased by a factor of four in the power 
spectrum. 
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FIG. 9. This figure illustrates the predicted uncertainties [see Eq. ( |50[ )] in the power spectrum from the Bright Red Galaxy 
(BRG) subsample of the Sloan Digital Sky Survey (SDSS). The solid curve is a CDM power spectrum with Q,h — 0.25. The 
bins in k and the errors are determined by Eq. (pd). Also shown by the dotted curves are two features in the power spectrum 
at the indicated values of fc* and for the same value oi A — 0.23 [see Eq. (M)]. 

The sold curve in Fig. is generated assuming a Harrison-Zel'dovich primordial spectrum and a CDM transfer 
function. The points are in wavenumber bins with errors taken from Eq. (pfl). The two dotted curves are power 
spectra with features of the form that would be generated by resonant particle production. In the form of Eq. (|4^), 
the two examples have k^, — lO^^/iMpc"^ and fc* = lO-^ft-Mpc"^ and both have A = 0.23 (the result of the model 
with A = 1, N = 1, and m = 2 x 10^^ GeV). Clearly either feature would be detected. 
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FIG. 10. The angular power spectrum calculated for a featureless Harrison-Zel'dovich CDM model (solid curve marked 
H-Z) and models with features produced by resonant particle production with A = 0.23 and the indicated values of k, [see Eq. 
(list)]. The cosmological parameters are identical for all models. The shaded area indicates the expected uncertainty in the Ci's 
from the MAP satellite. The models have all been normalized to the same value of C2. 

Features in the primordial power spectrum will also affect the CBR temperature fluctuations. The CBR temperature 
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fluctuations can be expanded in spherical harmonics, ST jT = J2iJ2m^i"iyim{d,4') (2 < Z < oo and —l<m< I). 
The anisotropies are described by the angular power spectrum, C; = (|a/mp), as a function of multipole number I. 

The amplitude of the perturbation spectrum at a given wavenumber k contributes to the angular power spectrum 
in a range of /, so the effect of a narrow feature is spread over several values of I. The angular power spectrum is 
easily calculated using CMBFAST |37|. The result of a sample model calculations is shown in Fig. M. The angular 
power spectrum resulting from a featureless power-law spectrum (Harrison-Zel'dovich) is shown by the solid curve. 
Also shown for comparison are three other models with features in the primordial power spectrum of the form in Eq. 
(H), with A = 0.23 and K = IQ-^hM-pc'^, IQ-'^hM-pc-^^ , and IQ-^hM-pc'^. 

The CBR is sensitive to the power spectrum in the interval from about 10~^/iMpc~^ to 10~^/iMpc~^. For fc» < 
10~'^/iMpc~^, a feature with A < 0.23 will be too small to be detected because cosmic variance will be large at 
small I. For fc» ^ 10~^/iMpc^^, the largest effect will be near the first Zel'dovich peak and could have an effect on 
determination of the geometry of the universe from the location of the peak. If k^, ~ 10~^/iMpc~^, the largest effect 
will be in the region of the secondary peaks. If fc* is much larger than 10~^ft,Mpc~^, the effect will be out of the 
range of CBR sensitivity. 

A meaningful detection limit is hard to quote; the limiting value of A depends on fc* (among other things). In 
Fig. O we indicate by the shaded regions the values of A that will differ more than a 2a from the best-fit featureless 
power-law spectrum expected to be determined by SDSS and MAP. 

The feature around fc* = 3 x 10~'^ft.Mpc~"'^ in the SDSS limit is due to the feature falling between the widely 
spaced bins in that region. In the figure we have assumed the normalization will be set by CBR determinations, i.e., 
the normalization is a free parameter for the MAP limits, but not for the SDSS limits. In principle, the region of 
sensitivity could extend to fc* greater than 0.1/iMpc^^ if the evolution of the power spectrum in the mildly nonlinear 
regime was completly understood. 

The form of the limit from MAP around fc* ^^ 3 x 10^^/iMpc^^ is due to the fact that the spectral feature in that 
region falls in the region of the Zel'dovich peaks. 

For 10^^ < fc* ^ 10^^/iMpc^^, the smallest A that has a 2a effect is about 6 x 10"^.^ This would correspond to 
the limit 



A < 0.2/iV2/5 



(51) 
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FIG. 11. The shaded region will produce a 2ct change in the matter power spectrum with the expected sensitivity expected 
for the SDSS-BRG subsample and the MAP CBR satellite mission. 



^Here we are not saying that A 
not be a good fit. 



6 X 10 ^ would be cleanly detectable, only that a featureless power-law spectrum would 
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The correct way to interpret this hmit is that if it is violated, the feature in the spectrum produced by resonant 
particle production has the potential of being detected in large-scale structure surveys. 

Again, any limit on A will depend on fc*. For fc* near 10~^/iMpc~^ and A as small as 0.01, the effect of resonant 
particle production will be noticeable. For larger A it is possible that the feature could be noticeable in both CBR 
and large-scale structure surveys. 

VI. CONCLUSIONS 

For the most part, analysis involving the power spectrum assumes that it is a featureless power law. Occasionally 
one includes smooth variation in the spectral index, but usually retaining the assumption that the power spectrum is 
featureless. 

In this paper we presented a reasonably straightforward scenario for producing a spike (or several spikes) in the 
power spectrum. While analysis assuming a featureless power spectrum is certainly a reasonable first step, it would 
be wise to keep in mind that when doing so one is making an approximation about the dynamics during inflation that 
may not be true. 

Let us now elaborate about the particle physics motivations of our work. When speaking about a 'model of 
inflation', the phrase is actually used by the community in two rather different ways. At the simplest level a 'model 
of inflation' is taken to mean a form for the potential, as a function of the fields giving a significant contribution to 
it. In single-field models this is just the inflaton field (f). If one knows V{(f)) and the field value at the end of inflation, 
this allows one to calculate the spectrum of density perturbations. 

At a deeper level, one thinks of a 'model of inflation' as something analogous to the Standard Model of particle 
interactions g^ . One imagines that nature has chosen some extension of the Standard Model, and that the scalar 
fields relevant for inflation are part of that model. In this sense a 'model of inflation' is more than merely a specification 
of the the potential of the relevant fields. It will provide answers to some fundamental questions such as whether 
the relevant fields and interactions have already been invoked in some other context and, if so, what are the mass 
parameters and the couphng constants in the sectors coupled to the inflaton field(s). 

Of course, it would have been wonderful if inflation already dropped out of the Standard Model of particle physics, 
but sadly that is not the case. On the contrary, it seems very likely that the physics responsible for an inflationary 
stage during the early evolution of the universe is due to some particle physics model whose typical energy is much 
larger than the weak scale and whose field values during inflation may be larger than the Planck mass. 

Well studied extensions of the Standard Model, such as supergravity and superstring theories, have the Planck 
mass as a fundamental scale. The exact motivations and goals of these theories beyond the Standard Model might 
be different, but for many applications to cosmology they have several common features. For instance, they contain 
in the spectrum particles with Planckian masses and are formulated in extra-dimensions. These extra-dimensions 
are compact and smaller than the three large spatial dimensions. It is therefore possible to dimensionally reduce the 
system and obtain and 'effective' 3-1-1 dimensional theory leaving behind a tower of Kaluza-Klein (KK) states (pyrgons 
p9|) whose mass is of the order of the inverse size of the extra dimension. Since the extra dimensions are expected to 
have a size characteristic of the Planck length, these KK states therefore have Planckian masses. Furthermore, these 
states are largely degenerate, the level of degeneracy depending upon the geometrical structure of the compact space. 
It is therefore quite natural to expect that a large number of fermions with the same mass, approximately mpj, couple 
to the sector responsible for the inflationary stage. In the simplest approach, this degeneracy factor may be accounted 
for by the parameter N we have used in the previous sections. Since N may be easily of the order of 100, a detectable 
signature in the spectrum may be present even in the case of small Yukawa couplings. Thus, the observation of a 
spike (or several spikes) in the power spectrum of the density perturbation may help us in understanding the features 
of the inflaton sector and thus testing Planckian physics. In the absence of any other experimental signature, this is 
certainly very intriguing. 
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